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Abstract 



In this paper a relativistic quantum theory introduced by de Groot and 
Ruijgrok in 1975, is used as a quark model in momentum space. The complete 
spectrum, with the exception of the selfconjugate light unflavoured mesons, is 
calculated. The potential used consists of an one-gluon-exchange (OGE) part 
and a confining part. For the confining part a relativistic generalization of 
the linear plus constant potential was used, which is well-defined in momen- 
tum space without introducing any singularities. For the OGE part several 
potentials were investigated. Retardations were included at all places. Using 
a fitting procedure involving 52 well-established mesons, best results were ob- 
tained for a potential consisting of a purely vector Richardson potential and 
a purely scalar confining potential. In this way a vector confining is entirely 
induced by the Richardson potential. Reasonable results were also obtained 
for a modified Richardson potential. Most meson masses, with the exception 
of the 7T 3 the K and the Kq, were found to be reasonably well described by 
the model. 

PACS number(s): 12.40.Qq, 03.65.Ge, ll.10.Qr 
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I. INTRODUCTION 



In principle the properties of mesons and baryons should be correctly described by Quan- 
tum chromodynamics (QCD). However, apart from some lattice gauge calculations, this is 
practically impossible at the moment. As a replacement simple quark models, in which 
hadrons are viewed as bound states of constituent quarks, are quite successful (for a review 
see Q). The simplest are the nonrelativistic (NR) ones. The potential used here, normally 
consists of a Coulomb term to account for the perturbative one gluon exchange (OGE), and 
a linear potential with possibly an additional constant to incorporate the nonperturbative 
confining. These models work very well for the heavier charmonia and bottonia. For the 
lighter mesons however, it is clear that relativistic corrections must be included. Roughly 
speaking, this can be achieved in two ways. The most direct way is to replace all NR ex- 
pressions by their relativistic counterparts. Spin dependences, like spin-spin, spin-orbit and 
tensor couplings are included by hand. The second way, which, from a theoretical point of 
view is more consistent, is to start from a framework that is manifest Lorentz covariant from 
the outset. The most natural representation for such a framework, like the Bethe-Salpeter 
equations, is momentum space. Traditionally, however, and also because of the belief that 
it would be impossible to describe a confining potential in momentum space, the equations 
are normally transformed to configuration space. Since a few years [0-0], however, it has 
been realized, that there is no obstacle to define a confining potential in momentum space, 
even in the relativistic case. Therefore a growing interest has arisen to study quark models 
directly in this more favourable representation. 

This will also be the subject of the present paper. The theory introduced by de Groot 
and Ruijgrok @-12] will be used as a model to calculate the masses of all known mesons, 
with the exception of the selfconjugate light unflavoured ones. This Lorentz covariant theory 
is defined via a natural generalization of the NR Lippmann-Schwinger equation and does 
not require further specification in the course of its solution. It does not start from the 
Bethe Salpeter equations. The main difference is that in the intermediate states all particles 
remain on their mass shell, and that total three velocity rather than four momentum is 
conserved. Negative energy states are not included. Retardation effects are incorporated in 
a simple and unambiguous way. This is to be contrasted with the Bethe-Salpeter equations, 
where different three-dimensional quasi-potential reductions lead to different retardations 
(see e.g. Sec. (2.3) of [|TJ). The theory has proven to give the correct fine structure formulas 
for the hydrogen atom and positronium ||12|| . In Sec. |J a brief summary of the theory will be 
given. In Sec. [TT1| a number of quark-antiquark potentials will be discussed. A modification 
of the Richardson potential, to acount for the OGE, as well as a relativistic generalization 
of the constant potential is defined. An important feature of the mesons consisting of light 
quarks is the appearance of linear Regge trajectories. Their origin in the light of the present 
theory is discussed in Sec. |TV| The numerical method used will be described in Sec. [V], and 
its results will be further discussed in Sec. IVl The paper ends with some conclusions. 
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II. FORMULATION OF THE THEORY 



In this section a summary of the theory, as introduced by de Groot and Ruijgrok 
with the modifications made in JO, will be given 

A. The general framework 

A state a of a quark (mass mi) and an antiquark (mass m 2 ), can be characterized by 
(piAi,p 2 A 2 ), where p\ and p 2 are the four-momenta of the quark and antiquark, and Ai and 
A2 are their helicities. Both particles are supposed to remain on their mass shell, which 



means that p° = y |pj| 2 + m| = E i: i = 1, 2. The theory is constructed in such a way that 
in the interaction the three-velocity 

Pi + P2 (n - v 

Pi+p 2 

is conserved. This means that the quark-antiquark potential Vp a for a transition from an 
initial state a = (pi\i,p 2 \ 2 ) to a final state (3 = contains only non-zero elements 

if v' = v. In the center of momentum system (cms) this velocity conservation coincides with 
three-momentum conservation, i.e., pi = — p 2 = p and p[ = — p' 2 = p'. In this frame 
therefore the potential can be written as 

Vp a = ^y iA ^ )AlA2 (p',p). 

In the NR case the momentum dependence of a central potential appears in the form |q| 2 = 
|p' — p| 2 . In the relativistic case this expression must be replaced by a covariant one. Here 
the usual replacement |q| 2 — > — q 2 cannot be used because, due to the lack of four-momentum 
conservation, the loss of momentum qi — Pi — p\ by the quark, will in general differ from 
the gain of momentum q 2 — p 2 — Pi by the antiquark. Instead the following obvious and 
symmetrical substitution is made 

|q| 2 -> -Qi ■ 12 = |p - p| 2 - T(p',p), (2.2) 

where the term r, defined by 

r(p',p) = (E 1 -E' 1 )(E' 2 -E 2 ), (2.3) 

is responsible for retardation effects. The theoretical justification for this replacement is 



twofold. In the first place it can be shown [12| to be consistent with velocity conservation. 



The second and more practical justification is, that in the case of the Coulomb potential, 
r automatically generates the correct form for the Breit interaction (see |12|). In the equal 



mass case r = — (E — E') 2 , which is exactly opposite to the retardation used by Gross and 
Milana j|] and Maung, Kahana and Norbury 0. The difference in sign will give the wrong 
sign for the Breit interaction, which in turn will effect the fine structure of positronium. In 
0] it was shown that Eq. ( |2.3| ) gives rise to the correct positronium fine structure formula. 

The relativistic wave equation from which the mass M of the meson is to be solved is in 
the cms given by (h = c = 1) 
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[E 1 + E 2 -M}q XlX2 (p) + 



(2.4) 



A 'l A 2 

where the wave function \E , a 1 a 2 (p) is normalized as 

'17111712 

I \^x 1 x 2 {P)\- 

XlX'2 



171x1712 

E[E' 2 



dp' = 0, 



E / I*Wp) 



E\E 2 



dp = 1 



(2.5) 



and V = 4i7iii7i2W . The quantity W is introduced for convenience, because it reduces in 
the NR limit to the NR potential. In this limit Eq. (|2.4j) reduces to the NR Schrodinger 
equation in momentum space. 



B. Decomposition 



The interaction W used in this paper can be decomposed into a vector part Vy and a 
scalar part Vg, which is in the cms given by 



W(p',p)=u K (p[)v X2 (p 2 ) [ 7 « • 7 (2)/i Mp', P) 
+ ^ 1 H^V s (p',p)} v^iQuxM. 
Here the Dirac spinors u and v for the quark resp. antiquark, are defined by 

P^ 



(2.6) 



u\{P) 
v\(p) 



N 
N 



1 

2Xb 



X(A, 



-2Xb 
1 



V 

(-i)er 2 X*(A, -), 



with iV = y (E + m) / (2m) , b = p/(E + m), and x(\p/p) the helicity spinor with helicity 
A. For the two-particle helicity states we use the conventions introduced by Jacob and Wick 
13| . Potential (|2.6j ) partially decouples with respect to the states \p; JM; X1X2), which are 



defined by Eq. (18) of |L3|], giving 



(p';J'M';X[)^\W\ P ; JM; AiA 2 ) 
5jj'W'(A' 1 A^|H/ J (p',p)|A 1 A 2 ). 



(2.7) 



Because of conservation of parity, W further decomposes into two 2x2 submatrices, each 
having a definite parity. The subspace spanned by 



1 1 ^ , 1 !\ 
n:n +|-o:-o 



1*2) 



1 

V2 



2 1 2 
1 1 



) + 



2' 2 
1 1 
2' 2 



(2.8) 



has parity (— 1) J+1 . It contains the triplet J = I ± 1 states. The complementary subspace, 
spanned by 
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|Sl> = 75 



1 1. 



1 1 



l*2> 



V2 



2' 2' 
1 1 



2' 2 



(2.9) 



2' 2 



) 



II) 

2'2 ; . 



has parity (— 1) J and contains the J = I singlet and triplet states. Only in the equal mass 
case this subspace further splits into two lxl subspaces. Let 



vr 3 j = p'p 



m\m2 



{n i \W J {p',p)\n j ), n = s,t, 



E[ E2E1E2 



(2.10) 



then eigenvalue Eq. ( |2.4| ) can be cast in the form (suppressing the quantumnumbers J, M 
and s or t) 



[ El + E 2 -M] fi(p) + ]T / VijtfMitfW = 0. 

■ J 



(2.11) 



In appendix |A| explicit formula's for V^ J = (Vy)™ J + (Vs)™ J are given. The reduced wave 
function / is normalized to 



E / \m\ 2 d P =i. 

; JO 



(2.12) 



III. THE QUARK- ANTIQU ARK POTENTIAL 

The quark- ant iquark potential must contain an one-gluon exchange (OGE) to account 
for the short range, and a confining part for the long range interaction. It is generally 
believed that Voce should have a vector Lorentz structure, while about the confining part 
Vcon there is no consensus. Some believe that it must have a purely scalar structure, while 
others admit a mixture between scalar and vector coupling. We will adopt this last point of 
view. The potential can therefore be written in the form (see Eq. (|2.6f) ) 



Vy — VoGE + eVcON, 

V s = (l- e)^coN, (3.1) 

where e is the scalar-vector mixing of the confining potential. 

For Vcon the relativistic generalization of the linear potential, as described in ]7| , plus a 
constant potential (to be defined below) was used. This generalization is defined in a formal 
way and does not introduce any singularities. For the OGE two different potentials were 
used: the Richardson potential Jl4| and a modified version of this potential, both containing 



a running coupling constant (see Sec. [Ill A|) . The Richardson potential contains a linear 



part by itself. Therefore in this case (from now on denoted by I), e = was chosen, so that 
the confining in the vector direction is completely determined by the Richardson potential. 
The modified Richardson potential has no linear part. Therefore in this case (denoted by 
II) a nonzero e was admitted. 
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In all these potentials the NR momentum transfer |q| 2 was replaced by —q\ ■ qi (see 
Eq. (|2.2j )) and in this way retardation effects were included everywhere. For notational 
convenience, the quantity 

Q = yj(- qi ■ q 2 ) (3.2) 
is introduced. In the NR limit it reduces to |q|. 



A. The one-gluon-exchange: 
running coupling constant 



The renormalization scheme of perturbative QCD says that, for large momentum trans- 
fer, the running coupling constant a s as it occurs in the one gluon exchange (the factor g 
arises from color averaging) 



OGE 



4a s (Q 2 



3 2vr 2 Q 2 



(3.3) 



is given by (see also Eq. (B.2) of JL5|) 



a s (Q 2 ) 



log(X n ) 
1 - b n — h O 



(3.4) 



Here 



12tt 
(33 - 2n) ' 



6(153 -19n) 
(33 -Inf ' 



X n = log [Q/A ( '° 



MS. 



and n is the number of quarks with a mass smaller than Q. The subscript MS denotes that 
the renormalization is performed according to the modified minimal subtraction scheme. The 



connection between the different A™'s is given by Eq. (B.4) of |15| . The typical momentum 



transfer within a meson is on the order of one GeV, so in this region n — 3. Therefore, 
Eq. ( p.4|) with n = 3 is in many cases used as an approximation for all large momentum 
transfers. In addition the 63 term is almost always neglected. But this term is not small 
at all: in the Q-region from 1 to 5 GeV its contribution is about 25%. Even for very high 
momentum transfers its contribution is substantial ~ 15% for Q = 50 GeV. However, it 
appears that when A^4r rather than A^- is used, a fairly good approximation of Eq. (B .41) 
for large Q is obtained by putting 

at, « Y~ (3 - 5) 

(see the curve "standard approximation" of Fig. [l|). For Q = 5 GeV the deviation from 
Eq. ( |3.4j ) is ~ 7%, and for Q ~ 50 GeV there is no detectable difference. Also for smaller 
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Q the agreement is better, but of course in this region the validity of Eq. (|3.4j ) is doubtful. 
Nevertheless we think that these considerations show that there is no theoretical necessity to 
stick to the value of a 3 in Eq. ( |3.5| ): a small deviation from it also results in a good running 
coupling constant for large Q. 

For small positive Q-values Eqs. ( |3.4| ) and ( |3.5| ) diverge. To remedy this, Richardson 
|T4| proposed a potential in which the divergence is shifted to the origin by making the 
replacement Q 2 —>■ Q 2 + A 2 in Eq. 



UQ 2 ) = - „ °° , 5* (3-6) 



2 7r 2 Q 2 log[l + g 
a A 2 a 



2vr 2 Q 4 4tt 2 Q 2 



+ . . . for Q -> 0. 



The color factor | is absorbed in a . In Fig. |T] the running coupling constant, defined via Eq. 
( [Of ) with Vcon = Vr, is compared to the QCD formula for ao = |ct3 = 167r/27 = 1.862. The 
singularity for Q = results from a linear term in the potential with string tension ^aoA 2 
(see or Sec. |III B|) . When the singularity is subtracted, the running coupling constant 



saturates to the value ^a^ = 0.698. From Fig. [I] it is seen that for momentum transfers 
starting from 2 GeV, a much better approximation to Eq. Q3.4Q is obtained, if ao is slightly 
decreased. A value of a = 1.750 turns out to be a very good choice. 

A different way to remove the singularities is to also make the replacement Q 2 —>■ Q 2 + A 2 
in 1/Q 2 itself. This results into 

V M {Q 2 ) = — nr- (3.7) 

27r 2 [Q 2 + A 2 ]log[l + f] 

n " +...forQ^0. 



27r 2 g 



This modified Richardson potential Vm does not contain a linear part. The coupling constant 
saturates to a value a 3 . The running coupling constant defined via this potential is given 
in Fig. [I] for ao = §03. From this figure it is seen that this choice for «o gives a good 
representation of the QCD formula for moderate Q values. 

The spinless partial waves Wr and Wm of Vr and Vm, defined by Eq. (|AT|), are given in 
appendix O. 



B. The confining: 
"linear + constant" potential 

For the confinement a relativistic generalization of a linear plus constant potential was 
used: 

Vcon = ^lin + Vc- 

As was already mentioned in the introduction, it was for a long time believed that a linear 
potential could not correctly be described in momentum space. A naive consideration shows 
that it behaves like —1/Q 4 , which results in an ill-defined bound state equation. A few 
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years ago it was shown 0,0,|5| that this singularity for the NR case is only apparent. For 
the relativistic case different methods were employed [f|,P,|7j to solve this problem. The one 
used in this paper is defined by [7| 



LIN 



T74.0 drf 2tt 2 



1 



Q 2 + V 2 



(3.8) 



where the color factor | is absorbed in the string tension A. In [7| it was shown that the limit 
exists in a distributional sense. The result was that the integral in Eq. ( 2.1 1|) is replaced by 



POO 

Jo 



vS iJ \p',p)W) + 



y 



4p 2 Ck(p) 
(p' 2 — p 2 ) 2 



dp 1 . 



(3.9) 



Here V- 



ii j 



n 



is the naive pointwise limit obtained from Eq. (|3.8| ). The l/(p' — p) 



singularity is removed by the quantity 



lim 



-\p'-pM nJ \p\p) 



(3.10) 



The resulting l/{p' — p) singularity is handled by the principal value integral (denoted by 
-p). It was shown that this subtraction is not just a trick to avoid singularities, but is really 
generated by Eq. ( |3.8| ). 

For a confining potential that consists of a mixture of a scalar and vector Lorentz struc- 
ture (see Eq. (|3.1|)), the pointwise limits of the spinless partial waves (see Eq. (|A1|)) are 



given by W v = eW[ and W l s 



(1 — e)W L , where 
W l L (p>, P ) = ±R(p',p)^ 

71 p'p 



(3.11) 



and also R(p',p) is given in Appendix Here z is defined by Eq. ( |B2|) and Qi is the 
Legendre function of the second kind of order /. The l/(p' — p) 2 singularity of W l L is 
determined by 



A 



R(p,p) 



The retardation defined by Eq. 

T(p',p) 



71 (p' — p) 2 — T(p',p) 
behaves around p' « p like 

2 



P 



(p'-p) +■■■ 



and therefore contributes to the singularity. It follows that is given by 



CM = - 

71 



€+(!-€ 



Pi - Pi 



(3.12) 



where p\ ■ p 2 = E 1 E 2 + |p| 2 is the dotproduct between the four vectors p\ and p 2 . Note 
that Cij does not depend on J and the parity s or t. In addition it is a manifest Lorentz 
covariant object. 
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When the interaction does not contain a linear part, the integral ( |3.9|) coincides with 
the integral within Eq. ( [2. 11 ). This is so because then = and there is no l/(p' — 
p) singularity, which means that the principal value coincides with an ordinary integral. 
Therefore replacement (|3.9| ) in combination with Eq. ( |3.10|) can be applied to the entire 
interaction. In this way a nonzero value of C automatically indicates the presence of a 
linear term. For the Richardson potential (|3.6|) with a purely vector character this results 



in 



RJij 



a A 2 
2vr ih 



(3.13) 



which indicates a linear term with string tension ^a^A 2 . 

In analogy with the linear potential the constant potential Vc can also be defined via 
the Yukawa potential. In the NR case one has in configuration space 



V c (r) 



C 



lim — — 

vio or] 



Ce 



;/r 



Therefore in momentum space an obvious relativistic generalization is 



V c (Q 2 ) = lim^- 
■nio or] 



C 



27r 2 (Q 2 + r] 2 ) 



Note that this expression also includes retardations, which are hidden in Q 2 (see Eq. 
Definition ( p,14|) has to be included in the integral of Eq. (|2.11|) before the limit is 
The spinless partial wave of this constant potential is given by 



(3.14) 

3))- 

taken. 



W l Jp\p) 



CR(p',p) r] 



71 



pp 



z + 



2p'p 



The only term that survives the limit r] j is 
CR 7] 



7T (p' — p) 2 — T + T] 2 



CR 



E\E 2 
Pi -P2 



6(p' -p) + - 



For a confining potential that has both a scalar and vector part (see Eq. ( |2.6|) ) this results 
in 



{V c )l J = C 



Pi -P2 



E 1 E 2 



e+(l-e 



mim2 I 6ij6(j/-p). 
Pi -P2 



(3.15) 



IV. LINEAR REGGE TRAJECTORIES 

The mesons which consist of the light u, d and s quarks only, are found to lie on so- 
called linear Regge trajectories. This means that there are several groups of mesons, for 
which the mass squared for each meson within such a group, is proportional to its angular 
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momentum J, i.e. Mj rs (3 J + C. The constant C depends on the group, the Regge slope 
(3 however is about the same for all groups. Its experimental value is (3 ~ 1.2 (GeV) 2 . For 
mesons containing a heavy c or b quark, such trajectories have not been observed. This 
makes it plausible that linear trajectories are induced by relativistic effects. In fact, it is 
known [|i^ |i~6| , |T7f , that the Schrodinger equation with ultrarelativistic (UR) kinematics (i.e. 
2p instead of p 2 /(2/z)) for a linear potential, does indeed give rise to linear trajectories, while 
the (NR) Schrodinger equation does not. The slope f3 solely depends on the string tension 
A, namely (3 — 8A. 

For the present case a similar effect is observed. It numerically appears that the (UR) 
limit (i.e. mi, m 2 —>■ 0) of bound state equation ( [2.1 1| ) also leads to linear trajectories, with a 
group independent slope (3. This slope however, depends on the vector part Ay of the string 
tension only. In addition, the dependence is a factor v2 larger than for the relativized 
Schrodinger equation, namely 



(3 « (8^2)A; 



(4.1) 



As can be deduced from Appendix [A], the off-diagonal elements Vu and V21 of both a vector 
and a scalar potential vanish in the (UR) limit. Therefore Eq.( [2.11| ) further decouples into 
two single equations. For the pure vector case, it reduces for the V// channel to 



[2p - M] fj(p) 



roo 
JO 



V J (p',p)fj(p') + -fj(p) 

7T 



dp' = 0, 



(4.2) 



with 



2A 



V J (p',p) = —Q[ 



Tip'P 



p + p + ip' — pY 

2p'p 



(4.3) 



This equation was solved numerically using the method described in Section 0. The calcu- 
lated masses (in units of vA) of the lowest states for each J are presented in Table |. The 
Schrodinger equation with (UR) dynamics is in momentum space also given by Eq. Q4.2|) , 
but now with 



A 



V j (p',p) = —Q: 



Tip P 



p' 2 + p 2 

2p'p 



(4.4) 



The corresponding calculated masses are also listed in Table |. They all agree with the 
calculations performed by Basdevant and Boukraa (see Table I of [|i~6|j). 

In principle the trajectories are expected to be linear only for large values of J. However, 
as Table | shows, the convergence is very fast. It was found that also for moderate masses 
Eq. (|2.11| ) leads to Regge trajectories with the same relation (|4.1|) between (3 and A. The 
convergence, however, is then slower. When in addition a OGE term and a constant are 
added to the potential, relation (|4.1| ) is affected. The change, however, is not very large. 



Therefore it can be concluded that, in order to obtain reasonable Regge slopes, the string 
tension in the vector direction Ay should be around Ay ~ 0.1 GeV 2 . 
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V. NUMERICAL METHOD 



The present model embraces eigenvalue equation ( |2.4j ) in combination with a quark- 
antiquark potential W consisting of an one-gluon-exchange part Voce and a confining part 
Vcon- The way in which Voce and Vqon enter in eigenvalue equation fl2.4| ) is given by Eqs. 
Q2.6|) and (|3.1| ). The OGE potential is determined by two parameters and A and the 
confining potential is determined by a string tension A and a constant C. Furthermore a 
parameter e can be introduced to give the confining potential a mixed scalar- vector character. 

The numerical solution of the model can be divided into two parts. The first one con- 
cerns the calculation of the masses of the mesons, from eigenvalue equation ( |2.4j ), given all 
parameters of the potential under consideration, and the quark masses. The second part is 
the fitting to the experimental data. The eigenvalue equation was solved by expanding the 
wave function into cubic Hermite splines (see [18[]). The integration region p e [0, oo) was 
projected onto the finite interval x G [0, 1] by x = (p — Po)/(p + Po), where p was chosen 
in the physical region. On this interval iV equidistant spline intervals were chosen on which 
2N spline functions were defined. The matrix elements of the resulting eigenvalue equation 
for the expansion coefficients only involved single integrations of the potential times a spline 
function. This is a major advantage of the spline method compared to the more conventional 
expansion techniques, where the evaluation of matrix elements involves two-dimensional in- 
tegrals. The integration was performed using Gauss-Legendre quadratures. In the case 
where the singular point p' = p was inside the region of integration, special care had to be 
taken. In the first place, an even number of abscissas centered around p' = p was used. In 
that way the principal value, which occurs for the confining potential, is automatically taken 
care of |5|jl9[|. Secondly, the logarithmic singularity ~ log(|p'— p\), which is induced by both 
the Coulomb and the confining potential, was separately handled by means of Gaussian 
quadratures based on a logarithmic weigth function (see e.g. Table 25.7 of p0|). Another 
important advantage of using Hermite splines, is their small nonzero domains. Therefore 
on each spline interval only a few of these splines (four for the Legendre and three for the 
logarithmic quadrature) were needed to obtain high accuracies. The matrix equation was 



solved using standard techniques pi| , giving the meson masses M;. 

The choice of the projection parameter po and the number of intervals N, depended on 
the specific meson. For instance, the typical momentum transfer for the T mesons (bb) is 
about 1 GeV, so po = 1 GeV. The masses of these mesons are all known to a high precision 
and they are all radial levels of the J PC = 1 channel. The T„ is the tenth radial state 
(n = 10). Therefore 20 spline intervals were needed to guarantee accurate results. The 
K\ however, is the only known J p = 2 + strange meson, apart from the unconfirmed K%'. 
Therefore N = 8 was sufficient to obtain reliable results, po = 0.5 GeV was a proper choice 
for this meson. 

The second part of the problem was to get a good fit to the experimental data. For this 
purpose the merit function 



X 2 (a 1 , .,a n ) = 



M^(a u .,a n )-Mr 12 



(5.1) 



has to be optimized with respect to the parameters ai, ., a n . Here i labels the mesons, Mf^ 
and M* he denote their experimental and calculated masses, and a, their weights. A nonlinear 
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regression method, based on the Levenberg-Marquardt algorithm was used to perform the 
fits (see section 15.5 of |21|). This method requires as input the explicit knowledge of 
the derivatives of the calculated masses with respect to the fitparameters. For the present 
complex situation this information is not known. It is only known that the derivatives of a 
meson mass with respect to quark masses it does not contain, is equal to zero. Therefore the 
derivatives where approximated in the least time consuming way by the following expression: 

dM^ M^(a u aj + A, ., a n ) - M^(a u a n ) 
~da~ A ' 

In this manner all required information, e.g. M* he and <9M* he /da>j, is obtained by calculating 
all meson masses (n + 1) times. A more sophisticated method would considerably increase 
this number. Approximation (|5.2|) turned out to be very effective: starting with a physically 
sensible set of parameters, after four or five steps convergence to an optimum was reached. 
The value of the parameter A appeared to be of minor importance. A = 0.04 (dimensionless, 
in GeV, or GeV 2 , depending on the dimension of aj) was found to be a good choice. 

All mesons regarded to be established by the 1992 Particle Data Group (in Table 
[TV] indicated by a "•") were used in the fit, with the exception of the selfconjugate (i.e. 
Isospin 0) light unflavoured ones. For a fair description of these mesons, an annihilation 
interaction from initial qq states to final q'q' states should be included. Also the charmed 
strange D* and D sJ were excluded, because of the uncertainty of their quantum numbers. 
Furthermore the up and down quark were considered to be of equal mass. In addition the 
electromagnetic interactions are completely neglected. Therefore the 7r° and n ± , the K° and 
K ± , and so on will be degenerate in this picture. Because of the indistinguishability of the 
u and d quark, from now on such a quark will be denoted by "v/d" . This accumulates to 
a total of 52 mesons: 11 light unflavoured (u/du/d), 11 strange (su/d), 4 charmed (cu/d), 2 
charmed strange (cs), 10 charmonia (cc), 2 bottom (bu/d) and 12 bottonia (bb). 

For the weight ai the maximum of the uncertainty dM° xp of the measured mass and 
the predictive power of the model, was taken. It is difficult to give an estimate for this 
predictive power. In the first place quark models have a phenomenological nature; there is 
no direct link with QCD. In the second place, the mesons are not stable particles, but in 
fact resonances. The decay mechanisms, which are not incorporated in this paper, could 
considerably effect the position of the calculated masses. This especially applies for the 
mesons that have decay widths of a few hundred MeV. To account for all of this, a grid size 
S = 20 MeV was introduced to give a minimum to a. Only for bottonium (bb), a grid size 
of 10 MeV was used, because in this system relativistic effects are less important and most 
states have narrow widths. Summarizing, the weights were determined by the following 
formula 

a { = Max [dMf p , S] . (5.3) 

A few exceptions to this rule were made. The pion 7r and the kaon K are the ground 
states of the u/du/d resp su/d mesons. It is commonly believed, that, in order to give a fair 
description of these particles the mechanism of chiral symmetry breaking should be included 
in the model. It appeared that also the Kq mass was badly described by the model. This 
state, however, has a large decay width of ~ 300 MeV. Therefore a n = 0.4 GeV and 
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&k = &k* = 0.2 GeV were chosen, so that these states get an insignificant weight in the fit. 
Another point are the p' and p" . These states also have large decay widths (~ 300 MeV). It 
appeared that best results were obtained if each state was regarded to be composed of two 
neighboring resonances. In Sec. [VI| his point will be discussed in more detail. 

To decrease the computation time first a rough fit was made by taking only half of the 
spline intervals N needed to obtain the desired accuracies. The resulting fitparameters were 
then used as the starting point for a full accuracy fit. The typical computation time for a 
complete rough fit for all 52 mesons was 30 min on a Sparc 2 workstation, while the fine 
tuning fit took about one hour. 

As was already mentioned, two different types of potentials were examined. In case I the 
Richardson potential Vr was taken to account for the OGE and for the confinement in the 
vector direction. Vcon has a purely scalar character (e = 0). For «o both the "QCD" value 
167r/27 and the value 1.75, which gives a better agreement with the QCD formula G3.4|) , 
were taken. Both choices ended in comparable fits (x 2 ~ 260). The resulting parameters 
for «o = 1-75 (denoted by la) are given in Table |T| and the calculated meson spectrum in 
Table [TV| Also the case in which naither ao nor e was fixed was regarded. The regression 
method led to very small values for e and the string tension A. Therefore a fit, denoted by 
lb, was made where these two parameters were put equal to zero, and where ao was varied. 
This resulted in a somewhat better fit (lb) with \ 2 = 250 (see also Tables [TTT] and |TV| ). In 
both cases seven parameters, three to model the potential and four quark masses, were fit. 
Finally, the case was considered in which the linear term of Vr was subtracted. To get a 
confining potential in the vector direction, the mixing e was also varied. The results did not 
improve, however. 

In case II the modified Richardson potential Vm in combination with a mixed scalar- 
vector Vcon wa s taken. The value «o = 167r/27 was the only parameter held fixed, so that 
eight parameters were varied. In spite of the extra parameter, the resulting fit (II), see 
Tables |ITJ and [TV], is worse than the fits found for case I and gave x 2 = 322. 

VI. DISCUSSION 

The meson spectrum calculated for parameter sets la, lb and II is given in Table [TV]. 
Also the mesons that were not involved in the fitting procedure (the ones without a a) 
were calculated. It is seen that most of these unconfirmed mesons (see |15[), are reasonably 
described by the model. Many states are a mixture between two 2s+1 Lj waves. Only in 
the NR limit these waves decouple because only then the angular momentum I is a good 
quantum number. For each state the most dominant wave is underlined. Most distributions 
are like 99% vs. 1%, which supports the statement that I is almost a good quantum number. 

A few years ago (for a review, see page VII of flT5l), the p(1450) and p(1700) were recog- 
nized as being a splitting in the formerly known p(1600) resonance. It could be interpretated 
as the fine structure splitting between the n = 2, dominantly S, and the n = 3, dominantly 
D, states. The splitting however is rather big (~ 250 MeV). For the present model this 
interpretation was found to be in conflict with the rest of the spectrum. A correct p' — p" 
splitting induced a far too large splitting in the 1 states of charmonium and bottonium, 
and visa versa. Only if the p' was regarded to consist of the n = 2 and n = 3 states, and 
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the p" of the n = 4 and n = 5 states, correct splittings for the entire spectrum could be 
obtained. In addition, the correct splitting between the observed p 1 " and p m (~ 50 MeV) 
was obtained. The difference between the n = 2 and n = 3 mass, and between the n = 4 
and n = 5 mass, was found to be ~ 100 Mev, which is much smaller than the decay widths 
of the p' and the p" (~ 300 MeV). 

The masses of the ir, K and Kq are found to be much to high. As was already mentioned, 
this was to be expected, because the small masses of these particles are believed to be a 
consequence of spontaneous breaking of chiral symmetry. In all cases la, lb and II, the mass 
of the rj c is found to be too high. As was pointed out by Hirono et.al. |22||, this is problably 



a consequence of neglecting negative energy states. They found that for a quark model 
based on the instantaneous ladder BS equation for charmonium, the r\ c is strongly influences 
by neglecting these states (~ 100 MeV), while the dependence on all other states is much 
weaker (~ 10 MeV). If one extrapolates these results to the present theory, this means that 
the omission of the negative energy states only weakly affects the spectrum. Only for the 
1 S'o ground states a substantial mass drop may arise. This would mean that also the masses 
of the D and the D$, which are now found a bit too high, would become smaller. The B 
would also get a smaller mass, which only has a positive result in case la. 
The centre of gravity COG(n) (see e.g. Sec. 8.1 of (TJ) is defined by: 

COG(n) = -M{n 3 P 2 ) + -M(n 3 Pi) + -M{n 3 P ) 
9 3 9 

It can be proved that, for an arbitrary scalar potential Vg and a Coulomb vector potential 
Vy, up to first order relativistic corrections, this COG equals the corresponding n^Pi singlet. 
The relation is violated by the Q-dependence of a s and the presence of a confining term in 
the vector direction. It is also affected by higher order relativistic corrections. In all cases 
the COG is found to be somewhat higher than the corresponding singlet state. A related 
quantity is the ratio p|.|24| 

9 M( 3 P0 - M( 3 P ) ' 

Its experimental value is 0.21 for u/du/d, 0.48 for cc, 0.66 for n = 1 bb and 0.57 for n = 2 
bb. For all three cases la, lb and II, a rather constant value of p ~ 0.8 (see Table [TIT]) was 
found. A perturbative configuration space calculation shows (see e.g. Sec. 4.2 of 0) that 
this too large value for p is a consequence of the dominance of the vector OGE. An analysis 
for the present case in momentum space, gives a similar result. A more profound linear 
scalar potential might lower this ratio. 



The following remarks on the parameter sets can be made. From Table |I| it is seen 
that the quark masses are substantially larger than usual in quark models. Furthermore, 
the masses are quite different for the different cases. The smallest masses are obtained by 
fit la. This is a consequence of the large negative constant C ~ —1.0 GeV, which, however, 
is neccesary in order to obtain a good fit for the entire spectrum. If, for instance one only 
considers bottonium and charmonium, it turns out that the quality of the fit only weakly 
depends on the value of C. The system is overparametrized and one in fact does not even 
need a constant in the potential. But, when simultaneously also good results for the lighter 
mesons are required, the large negative constant arises automatically. 
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The total string tension A to t is defined as the sum of the tensions in vector and scalar 
direction. For case I one has A to t = A + |a A 2 , while case II simply gives A to t = A. These 
tensions are also quite different for the different cases. In case lb, which gave the best fit, 
there is only a vector tension. This is in contrast to the requirement that, in order toobtain 
better p values, the confining should be dominantly scalar. The total tension for case la is 
closest to value A ~ 0.18 GeV 2 , which is often given in the literature. 

The Regge slopes of la and lb are compatible with the experimental value f3 ~ 1.2 GeV 2 . 
The slope found in II, is somewhat too low. This can clearly be seen from the high- J states 
like the ps, a 6 and K§. The errors given in Table [Til] represent a measure of linearity of 
the trajectories. It is defined as the spread in the difference between the masses squared of 
adjacent states. The spreads found are considerably smaller than the experimental value. 

Finally the running coupling constant for Q = 31 GeV and Q = Mz = 91.16 GeV were 
compared with the experimental values 

ct s (34 GeV)= 0.14 ±0.02, 
a s (M z ) = 0.1134 ±0.0035. 

Only case la is compatible with both conditions. The choice ao = 1.75 was made to give 



the best approximation to the QCD formula ( |3.4| ) for moderate momentum transfer. Now 



it appears that this choice also gives correct results for very high momenta. A fit of type la, 
but now with a® = 167r/27 (not displayed) gave a too large high momentum a s . In principle, 
the high Q-range of the potential is completely irrelevant for the calculation of the meson 
spectrum, where the potential is only tested up to a few GeV. Nevertheless, for the sake of 
theoretical consistency, this test was made. 



VII. CONCLUDING REMARKS 

In this paper a relativistic quark model defined in momentum space was studied. The 
quark- ant iquark potential used, consisted of a OGE with a Lorentz vector character, and a 
linear plus constant confining potential. For the OGE the Richardson potential Vr, given 
by Eq. (|3.6|) , with and without its linear part, as well a modified Richardson potential Vm, 
defined by Eq. (|3TT|), was regarded. Best results were obtained for the Richardson potential 
including its linear term (case I). The linear plus constant potential was given a pure scalar 
character, i.e. e = in Eq. ( |2.6j ). In this way, the confining in the vector direction was 
induced by the linear part of Vr. For case I two different fits were made, fit la, in which 
the value of «o was fixed to 1.75, and fit lb, in which «o was varied, but the string tension 
in the scalar direction A was put equal to zero. Also reasonable results were obtained for 
Voge = Vr (case II). Here the confining potential was given a mixed scalar- vector character. 
For the fits la, lb and II, most meson masses, with the exception of the tt, K and Kq were 
found to be reasonably described by the model. In case la and lb correct Regge slopes were 
found, and only in case la a correct strong coupling constant for large momenta was found. 
The ratios p, defined by Eq. (|6.1| ), however, were in all three cases found to be too large. It 
is concluded that case la should be preferred. 

No detailed comparison with other theories has been made, because the main purpose 
of this paper was not so much to improve upon the existing calculations, but rather to 
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show that results of the same quality could be obtained using a relativistic theory which is 
formulated in the momentum representation. 
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APPENDIX A: PARTIAL WAVE DECOMPOSITION 

In this appendix we give the precise form of the decomposition of a potential W defined 
by Eq. (|2.6|). The partial waves, 

= (Vv)# + (V 3 )lf, n=s,t, 

defined by Eq. ( |2.10|) can in general be expressed in terms of the "spinless" partial waves 
Wy and Wg of Vy and Vg respectively. They are defined by 

KAP'>P) = (Wp)R(j/,p) J^W^sip'^P^dx, (Al) 

with x = and Pi the Legendre polynomial of order I. The quantity R is defined by 
R(p',p) = A l N 2 A x A 2 , with A = J^tp. If furthermore b = -g^-, then the result is: 

1. Vectorpotential for the |si,2) states 

(V V ){{(p', P ) = [1 + 3(b[b' 2 + b,b 2 ) + b'^bM + (b[ ~ b' 2 )(h - b 2 ) { ^ ' 2 V J+ [ V , 

i W J+1 — W J ~ l 

{Vv)\l(p',p) = (b' 2 - b[)(b 2 + h^JiJ + l) \ J + 1 V = (V v y 2 ((p,P), (A2) 

TW J+1 + ('7 4- DW -7-1 

iv v y 2 i{p',p) = (i + b[b' 2 )(i + hb 2 )w v + (b[ + b> 2 )(h + b 2 ) v + [ ) v . 



2. Vectorpotential for the |ii 2) states 



(v v y 1 J 1 (p',p) = (i-b[b' 2 )(i-b 1 b 2 ) 



2J+1 

(Vv)%(jt,P) = "(I - b[b' 2 )(l + ' = (VvY 2 i(p,p'), (A3) 



+ [(b'i ~ 6s) (61 - 62) + 4(b[b 2 + b' 2 b x )\ W( 



TW J+1 4- ( T 4- 1 "lW/ J ~ 1 
(^)M(P',P) = (1 + &'l& 2 )(l + 6162) V or ,T + ( & i + + h ^ W V- 



2J + 1 

] 

2J4l" 
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3. Scalarpotential for the \si^) states 

(j + \)wi +1 + jwi- 1 



(Vs) s A(p',p) = (1 + h'&hb^W's - {b% + 6 2 6 2 )- , } ; i 

/ w J+1 — w 1 ^ 1 

(Vs)tiip',p) = (Kh - b' 2 b 2 )^J{J+l) S 2J + 1 S = (Vs) s 2 i(p, P '), (A4) 

TW J+1 4- ( T 4- 1 "llV J_1 

0^(p',p) = (1 + WiM^/ - (Kh + 5 2J + 1 - 



4. Scalarpotential for the |ii )2 ) states 
(V s )[i(p',p) = (1 + &;& 2 M 2 ) (J+1)1 !f'| J ^" 1 - (6i6i + b' 2 b 2 )W J s 



2J+ 1 

/ W J+1 — w J ~ x 

(V s )\i(p',p) = -(1 - MWy/JV+l) S 2J+1 S = <Ys)%M> (A5) 

(V^OAp) = (1 + WiM JW " +1 + 2 j J + \ l)WJs ~ l - ^ + lfa)Wi. 

Strictly speaking, these results are only valid for J > 0. For J = only the Vn's are 
nonzero and are also given by Eqs. (|A^,.., [A5l) , but with W 1 ^ 1 = 0. 

In the equal mass case there is no difference between the bi's and the b 2 's. From this it 
is seen that the Vg's and the V 2 \ J, s are zero. This means that the potential decouples with 
regard to the \si) state, which corresponds to the / = J singlet, and the \s 2 ) state, which 
corresponds to the I = J triplet. Therefore in this case only the \l = J ± 1) triplet states 
mix. In the unequal mass case the I = J singlet and triplet states mix. 



APPENDIX B: SPINLESS DECOMPOSITION OF THE RICHARDSON AND 
MODIFIED RICHARDSON POTENTIAL 

In this appendix the spinless partial wave decomposition Wr and Wm of the potentials 
Vr and Vm-i defined by Eqs. ( |3.6| ) and (|3.7|) will be calculated. 

The momentum transfer Q 2 which is defined by Eq. (|3.2|) depends on the lengths p — Ipl 



and p' = |p'| of the incoming resp. outgoing momentum, and on the angle x = between 
these two momenta: 

Q 2 = Q 2 (p,p',x) = 2pp'(z - x). (Bl) 

Here 

p 2 +p' 2 - r(p,p') 
*b(p,p) = ^ , (B2) 

where the retardation r is a theory dependent quantity which in the present case is given 
by Eq. (|2.3|) The spinless partial wave W l (p',p) corresponding to an angular momentum / 
is defined by Eq. (|A1[ ). Introducing 
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V(x) = 1 + b = (B3) 



A 2 ' 2pp' 
and y± = ?/(x = =Fl)>l, then W l R and W 7 "^ are given by 



2tt A_ (t/-l)lo g! / ^ ' 

2tt V y logy 
The y dependence in P z can be expanded, using 

P l {z-w) = Y.g\{z)w\ (B6) 

Here g\ is a polynomial of degree I — %. For / < 3 it is given in Table I. For general I it can 
be found from the recurrence relation 

(I + l)gl +1 = (21 + l){zg\ - g\_ x ) - lg l r\ (B7) 

in combination with the initial values 

g \z) = 0, g° (z) = 1. (B8) 

Note that g l Q obeys the recurrence relation of the Legendre Polynomials. In combination 
with the initial values ( |B~8| ) it follows that g l Q = Pi. 

The partial waves W l R and W l M can be written in terms of s and the integrals 

K = -\ ^y^-dy, n > -1. (B9) 
2 J y _ ylogy 



For n = — 1 one has 



where 



A-x = ^ [F(logy+) - F(logy_)] , (BIO) 



roc df; 

F[X ^~L u7^Ty x>0 - (B11) 



For n > the integrals A n can be expanded into 



b 



n n 



where 



An=jj:^l)(-l) (n - k) h, (B12) 



4 = / f dy, n > 0. (B13) 

logy 
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One finds 



log 



log y- 



I n = Ei(n\ogy + ) - Ei(n log yJ), n > 0, 



where 



Ei(x) 



J — CO t 



dt 



(B14) 
(B15) 

(B16) 



is the Exponential integral (see Eq. (5.1.2) of [0). The principal value integral is denoted 

Dummanzme; all steps it follows that W l R and W l M , defined by Eqs. (^§ and (§77p 

are given by 



71 



i=0 



Wm(p.p') 



7T 



(B17) 
(B18) 



i=0 



where the polynomials g\ are defined by Eq. ( |B6|) and the integrals A n can be found from 
Eqs. (IBIOD and (|Blp . 
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FIGURES 

FIG. 1. Running coupling constant a s (Q 2 ), defined by Eq. ( |3.3| ) for three different 
choices of Voce compared to the QCD formula ( fDj ) and its standard approximation (|3.5|) , 
A = A§L = 0.3 GeV and a = 16vr/27. 
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TABLES 



TABLE I. Regge trajectories calculated from the ultrarelativistic Eq. (4.2). The masses are 
expressed in terms of vA, where A is the string tension. 







Present case: 


Basdevant and Boukra: 






potential (|4l 


1) 




potential (JO|) 


1 


Mi 


{Mf- 


Mf_ x )/(8V2) 


Mi 


{Mf - M^/8 





3.830 






3.157 




1 


5.062 




0.969 


4.225 


0.985 


2 


6.066 




0.987 


5.079 


0.994 


3 


6.931 




0.993 


5.811 


0.996 


4 


7.701 




0.996 


6.461 


0.998 


5 


8.402 




0.998 


7.052 


0.999 


6 


9.049 




0.998 


7.597 


0.998 



TABLE II. The polynomials g\ for / < 3, defined by Eq. ( |B6j ) 



1 

1 z -1 

2 2" Z ^ ^Z 2" 

o 5„3 _ 3„ _15~2 , 3 15„ _5 

2 2 2^~ r 2 2 ^ 2 
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TABLE III. Final parameter sets from the fitting procedure described in Sec. |V] for potential 
models I and II. The varied parameters are indicated by a "•" . For model I two different fits were 
made. In case la ao was held fixed and A was fitted, while in case lb A was put equal to and ao 
was fitted. The related quantities are discussed in Sec. [V]. 



Model: 


la 


lb 


II 


Potential 








a 


1.750 


2.434 • 


1.862 


A (GeV) 


0.324 • 


0.320 • 


0.376 • 


A (GeV 2 ) 


0.077 • 





0.136 • 


C (GeV) 


-1.297 • 


-1.291 • 


-1.038 • 


e 








0.523 • 


Quark masses 








m u / d (GeV) 


0.512 • 


0.699 • 


0.966 • 


m s (GeV) 


0.766 • 


0.889 • 


1.072 • 


m c (GeV) 


2.066 • 


2.206 • 


2.249 • 


m b (GeV) 


5.474 • 


5.616 • 


5.593 • 


# parameters 


7 


7 


8 


x 2 


263 


250 


322 


Related quantities 








A tot (GeV 2 ) 


0.169 


0.125 


1.136 


(3 (GeV 2 ) 


1.18±0.05 


1.27±0.09 


0.93±0.08 


P 


0.81 


0.79 


0.75 


a s (34GeV) 


0.141 


0.196 


0.155 


a s (M z ) 


0.1164 


0.161 


0.127 
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TABLE IV. Meson spectrum calculated from Eq. (|2,4j) for three different parameter sets la, lb 
and II, ( see Table pi] ). All masses are in (GeV). The experimental values are taken from [15|, with 
the exception of the h c \, which is taken from p3| . The mesons labeled with a "•" (regarded as 
beeing established by fl5|]) were, with the exclusion of the D* and the D s j, involved in the fitting 



procedure. 


The weij 


jhts <Jj are 


determined by Eq. (5.3) 


. The most dominant 


2s+l Lj 


waves are 


underlined. 


















Light unflavoured mesons: u/d quarks. 


Name i 


J PC 


2s+ i Lj 


M exp 


n 


my 


Mf 


Mf 




• 7T 


o-+ 


1 5'o 


0.135 


1 


0.600 


0.595 


0.688 


0.400 


• 7T' 


o-+ 


1 5'o 


1.300 


2 


1.243 


1.206 


1.292 


0.100 


7r" 


o-+ 


1 So 


1.775 


3 


1.711 


1.671 


1.695 




• P 


1— 


3 S 1 / 3 D 1 


0.768 


1 


0.754 


0.762 


0.867 


0.020 




1— 


3 S 1 / 3 D 1 




2 


1.365 


1.345 


1.387 




• p' 


1— 


3 5i/ 3 ^ 


1.465 


3 


1.474 


1.477 


1.460 


0.025 


• P" 


1— 




1.700 


4 


1.806 


1.786 


1.764 


0.020 




1— 


3 S 1 / 3 D 1 




5 


1.865 


1.864 


1.807 




P'" 


1~ 


3 S 1 / 3 D 1 


2.100 


6 


2.162 


2.151 


2.065 




p iv 


1— 


3 S 1 / 3 D l 


2.150 


7 


2.200 


2.206 


2.096 




• a Q 


0++ 


3 p 

-TO 


0.983 


1 


1.012 


0.981 


1.017 


0.020 




0++ 




1.320 


2 


1.517 


1.464 


1.510 






1++ 




1.260 


1 


1.166 


1.163 


1.197 


0.030 


• a 2 


2++ 




1.318 


1 


1.301 


1.319 


1.329 


0.020 


COG 




3 p 
0,1,2 


1.262 


1 


1.224 


1.229 


1.250 




• h 


1+- 


X Pl 


1.232 


1 


1.183 


1.194 


1.231 


0.020 




2-+ 


^2 


1.670 


1 


1.590 


1.614 


1.561 


0.020 


A 


2-+ 


^2 


2.100 


2 


1.958 


1.986 


1.880 






3~ 


3 ^3/ 3 C3 


1.691 


1 


1.698 


1.734 


1.637 


0.020 




3~ 


3 D 3 / 3 G 3 


2.250 


2 


2.051 


2.092 


1.940 






3~ 


3 D,/ 3 G, 




3 


2.097 


2.152 


1.969 






3++ 


3 F 3 


2.050 


1 


1.915 


1.957 


1.813 




a i 


4++ 


3 FJ 3 H A 


2.040 


1 


2.021 


2.085 


1.883 




P5 


5~ 


3 G 5 / 3 h 


2.350 


1 


2.297 


2.395 


2.093 




«6 


6++ 


3 K & / 3 .h 


2.450 


1 


2.540 


2.677 


2.279 





Strange mesons (Kaons): s, u/d quarks. 



Name i 


J PC 


2s + l Lj 


M exp 


n 


MY 


M} b 


Mf 




• K 

K' 
K" 


0" 
0" 
0" 


'So 


0.495 
1.460 
1.830 


1 

2 
3 


0.762 
1.402 
1.864 


0.723 
1.329 
1.786 


0.781 
1.385 
1.786 


0.200 


• K* 

• K*' 

• K*" 


1- 
1- 
1- 
1- 


3 S l / 3 D 1 
3 S 1 / 3 D 1 
3 Si/ 3 D 1 
3 S l / 3 D 1 


0.894 
1.412 
1.714 


1 

2 
3 

4 


0.891 
1.504 
1.618 
1.945 


0.876 
1.455 
1.588 
1.890 


0.955 
1.477 
1.553 
1.852 


0.020 
0.020 
0.020 


• K* 


0+ 


3 Po 


1.429 


1 


1.177 


1.112 


1.115 


0.200 
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0+ 


3 Po 




2 


1.674 


1.587 


1.604 






K*' 


0+ 


3 Po 


1.950 


3 


2.074 


1.989 


1.955 




• 


K x 


1+ 


'PJ 3 P1 


1.270 


1 


1.304 


1.274 


1.288 


0.020 


• 


K[ 


1+ 


'Pl/ 3 E1 


1.402 


2 


1.322 


1.306 


1.321 


0.020 




K'{ 


1+ 


'Ei/ 3 Pi 


1.650 


3 


1.773 


1.725 


1.701 




• 


K* 2 


2+ 


3 E 2 / 3 F 2 


1.429 


1 


1.416 


1.415 


1.415 


0.020 




K* 2 ' 


2+ 


3 E 2 / 3 F 2 


1.980 


2 


1.867 


1.849 


1.785 








2+ 


3 Pil 3 E_ 2 




3 


1.945 


1.934 


1.831 






K 2 


2- 


'D 2 / 3 D 2 


1.580 


1 


1.706 


1.693 


1.640 




• 


K' 2 


2~ 


'D 2 / 3 D 2 


1.768 


2 


1.715 


1.711 


1.650 


0.020 






2~ 


'D 2 / 3 D 2 


2.250 


3 


2.082 


2.065 


1.962 




• 




3~ 


3 D 3 / 3 G 3 


1.770 


1 


1.801 


1.813 


1.722 


0.020 






3+ 




2.320 


1 


2.027 


2.034 


1.900 




• 


K% 


4+ 


3 FU/ 3 H 4 


2.045 


1 


2.115 


2.150 


1.967 


0.020 




K 4 


4 _ 


1 G r 4/ 3 G 4 


2.500 


1 


2.300 


2.333 


2.116 






Kt 


5" 


3 G 5 /% 


2.380 


1 


2.386 


2.449 


2.176 





Charmed mesons: c, u/d quarks. 





Name i 


J PC 


2s+l Lj 


M exp 


n 




M} h 


Mf 


<7i 


• 


D 


0" 


'So 


1.867 


1 


1.935 


1.901 


1.904 


0.020 


• 


D* 


1- 


S £i/ 3 A 


2.010 


1 


2.006 


1.999 


2.031 


0.020 


• 


F>x 


1+ 


'EJ 3 Pi 


2.424 


1 


2.406 


2.379 


2.382 


0.020 




Dm 


1+ 


l Pi/ 3 P, 


2.440 


2 


2.439 


2.438 


2.424 




• 


D* 2 


2+ 


3 E 2 / 3 F 2 


2.459 


1 


2.485 


2.492 


2.484 


0.020 



Charmed strange mesons: c, s quarks. 





Name i 


jPC 


2s + l Lj 


M ex P 


n 


M} & 


Mf 


Mf 


0"i 


• 


D s 


O - 


'So 


1.969 


1 


2.032 


1.990 


1.984 


0.020 


• 


D*(?) 


1- 


3 S l / 3 D l 


2.110 


1 


2.100 


2.088 


2.110 




• 


D s i 


1+ 


'Ei/ 3 Pi 


2.537 


1 


2.498 


2.473 


2.466 


0.020 






1+ 


'Pl/'E! 




2 


2.520 


2.516 


2.503 




• 


D s j(?) 


2+ 


3 P 2 / 3 F 2 


2.564 


1 


2.561 


2.568 


2.563 





Bottom mesons: b, u/d quarks. 



Name i 




2s + l Lj 


M exp 


n 


Mf 


Mf 


Mf 


<?i 


• B 


o- 


'So 


5.279 


1 


5.303 


5.268 


5.247 


0.020 


• B* 


l - 




5.325 


1 


5.336 


5.318 


5.316 


0.020 



Charmonium: c quarks. 



Name i 


jPC 


2s + l Lj 


M exp 


n 


Mf 


Mf 


Mf 


0z 


• r] c 

rfc 


o-+ 
o-+ 


'So 
'So 


2.979 
3.590 


1 

2 


3.042 
3.615 


3.010 
3.589 


3.007 
3.609 


0.020 


• J/* 


1— 
1— 


3 £i/ 3 £>i 

3 S 1 / 3 D l 


3.097 
3.686 


1 

2 


3.099 
3.655 


3.104 
3.646 


3.117 
3.665 


0.020 
0.020 
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• 




1— 


s Si/ 3 £i 


3.770 


3 


3.766 


3.780 


3.775 


0.020 


• 


\$>"> 


1— 


3 S 1 / 3 D 1 


4.040 


4 


4.051 


4.017 


4.028 


0.020 


• 




1— 


3 V 3 £i 


4.159 


5 


4.124 


4.105 


4.097 


0.020 


• 




1— 


3 £i/ 3 A 


3.415 


6 


4.376 


4.319 


4.314 


0.020 






1— 


s Si/ 3 £i 




7 


4.430 


4.384 


4.364 




• 


XcO 


0++ 


3 J Po 


3.415 


1 


3.437 


3.433 


3.409 


0.020 


• 


Xcl 


1++ 


3 Pi 


3.511 


1 


3.485 


3.506 


3.504 


0.020 


• 


Xc2 


2++ 


3 £ 2 / 3 ^ 2 


3.556 


1 


3.523 


3.562 


3.572 


0.020 




COG 




3 p 
^0,1,2 


3.525 


1 


3.501 


3.529 


3.531 






Ki 


1+- 


'Pi 


3.526 


1 


3.492 


3.520 


3.522 




Bottonium: b quarks. 




Name i 


J PC 


2s + l Lj 


M exp 
i 


n 


M* a 


Mf 


Mf 


<?i 


• 


T 


1— 




9.460 


1 


9.493 


9.434 


9.441 


0.010 


• 


I 


1 


3 c /3 


10.023 


2 


10.011 


10.018 


1 n nil 
10.022 


0.010 






1— 


3 S 1 / 3 D 1 




3 


10.131 


10.171 


10.160 




• 


Y" 


\— 


3 £i/ 3 A 


10.355 


4 


10.346 


10.348 


10.365 


0.010 






1— 


3 5 1 / 3 a 




5 


10.423 


10.444 


10.451 




• 


Y'" 


1— 


3 £i/ 3 A 


10.580 


6 


10.614 


10.599 


10.626 


0.010 






\— 


3 V 3 £i 




7 


10.672 


10.670 


10.688 




• 


Y*i> 


\— 




10.865 


8 


10.846 


10.811 


10.844 


0.010 






\— 


3 5 1 / 3 D 1 




9 


10.893 


10.868 


10.892 




• 


Y^ 


\— 


3 £i/ 3 A 


11.019 


10 


11.054 


11.000 


11.035 


0.010 


• 


Xw 


0++ 


3 p 


9.860 


1 


9.859 


9.863 


9.843 


0.010 


• 


Xbo 


0++ 


3 J Po 


10.232 


2 


10.220 


10.229 


10.232 


0.010 


• 


Xbi 


1++ 


3 Pl 


9.892 


1 


9.882 


9.906 


9.888 


0.010 


• 


Xbi 


1++ 


3 Pl 


10.255 


2 


10.239 


10.258 


10.261 


0.010 


• 


Xb2 


2++ 


3 £ 2 / 3 ^2 


9.913 


1 


9.901 


9.938 


9.922 


0.010 


• 


Xbi 


2++ 


3 £ 2 / 3 ^2 


10.268 


2 


10.253 


10.281 


10.284 


0.010 




COG 




3 p 
^0,1,2 


9.900 


1 


9.890 


9.919 


9.902 






COG 




^0,1,2 


10.260 


2 


10.245 


10.267 


10.271 
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